Two-loop scalar self-energies and pole masses in a general 
renormalizable theory with massless gauge bosons 
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I present the two-loop self-energy functions for scalar bosons in a general renormalizable theory, 
within the approximation that vector bosons are treated as massless or equivalently that gauge 
symmetries are unbroken. This enables the computation of the two-loop physical pole masses of 
scalar particles in that approximation. The calculations are done simultaneously in the mass- 
independent MS, DR, and DR renormalization schemes, and with arbitrary covariant gauge fixing. 
As an example, I present the two-loop SUSYQCD corrections to squark masses, which can increase 
the known one- loop results by of order one percent. More generally, it is now straightforward to 
implement all two-loop sfermion pole mass computations in supersymmetry using the results given 
here, neglecting only the electroweak vector boson masses compared to the superpartner masses in 
the two-loop parts. 
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I. INTRODUCTION 

Low-energy supersymmetry provides a way of under- 
standing the small ratio of the electroweak breaking scale 
in the Standard Model to other very high energy scales. 
This requires the existence of complex scalar superpart- 
ners for each of the known quarks and leptons, with 
masses not far above 1 TeV. These squarks and sleptons 
should be accessible to the Fermilab Tevatron pp collider 
or the CERN Large Hadron Collider, and their masses 
can be measured with refined precision at a future linear 
e + e~ collider; see for example [1] and references therein. 
The match between these measurements and particular 
models of supersymmetry breaking will then require cal- 
culations at the two-loop level of precision, at least. An 
important part of this is the calculation of self-energy 



functions, which can in turn be used to calculate the 
physical masses of the new particles. 

In general, the mass given by the position of the com- 
plex pole in the propagator is a gauge-invariant and 
renormalization scale- invariant quantity [2]- [9]. The pole 
mass docs suffer from ambiguities [10] due to infrared 
physics associated with the QCD confinement scale, but 
these are probably not large enough to cause a practi- 
cal problem for strongly-interacting superpartners. The 
pole mass should be closely related in a calculable way to 
the kinematic observable mass reported by experiments 
[11]. In recent years, many important higher-order cal- 
culations of self-energy functions and pole masses in the 
Standard Model have been performed, including two-loop 
[12]-[15] and three-loop [16]-[17] contributions for quarks, 
and two-loop results for electroweak vector bosons [18]- 
[24] , as well as two- loop results for top and bottom quarks 
in supersymmetry [25]. 

In a previous paper [26] , I provided partial results for 
the two-loop self-energy functions for scalars in a gen- 
eral renormalizable theory, using the approximation that 
no more than one vector boson propagator is included. 
That is a useful approximation for the Higgs scalar bo- 
son(s), for which the most important contributions at 
the two- loop level involve the strong interactions and/or 
Yukawa couplings. In this paper, I will extend the previ- 
ous result by including the contributions for any number 
of vector boson lines, within the approximation that the 
gauge symmetry is unbroken so that the vector bosons 
are massless. Because of the experimental exclusions 
of light sfermions already achieved by the CERN LEP 
e + e~ collider [27] and the Fermilab Tevatron pp collider 
[28, 29], this will likely give a very good approximation 
for the squark and slepton pole masses. (Here, the ef- 
fects of non-zero W and Z boson masses can be included 
as usual in the one-loop part [30], and neglected in the 
two-loop part.) 
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II. NOTATIONS AND SETUP 

Let us write the tree-level squared-mass eigenstate 
fields of the theory as^ a set of real scalars two- 
component Weyl fermions ipj, and vector bosons V^. 
Scalar field indices are i,j,k,.. ., fermion flavor indices 
arc /, J, K, . . ., and A, B,C, . . . run over the adjoint rep- 
resentation of the gauge group, while /i, v, . . . are space- 
time vector indices. Repeated indices of all types are 
summed over unless otherwise noted. The kinetic part of 
the Lagrangian is taken to be: 

Adn = -^Rid^Ri-^mjR^ 

-i^a^d^j - ^(M'-V/VV + cc.) 

-\(d»Vt-d v V£)d*VZ. (2.1) 

The metric tensor has signature ( — h++). The non- 
gauge interactions of these fields are given by: 

Ant = —-z^ k RiRjRk — 7r7^ 4:,fem RiRjRkRm 
6 24 

~{y JKi MKRi + c.c). (2.2) 

where A'- 7 '* and A yfem are real couplings and the Yukawa 
couplings y JKl are symmetric complex matrices on the 
indices J, K, for each i. Raising or lowering of fermion 
indices implies complex conjugation, so 

M u = (M IJ y, y JKl = (y JK r. (2.3) 

The heights of real scalar and vector indices have no sig- 
nificance, and are chosen for convenience. The scalar 
squared masses m 2 and the fermion squared masses 
MjkM kj = ntfSj are taken to have been diagonalized 
(by an appropriate rotation of the fields if necessary). 
However, the fermion mass matrix M IJ is not necessar- 
ily diagonal, but must have non-zero entries only when 
I and J label two-component fermions with the same 
squared mass and in conjugate representations of the 
gauge group. 

In order to completely specify the pertinent features 
of the gauge interactions of the theory, let T A be the 
Hcrmitian generator matrices of the gauge group for 
a (possibly reducible) representation R. They are la- 
beled by an adjoint representation index A correspond- 
ing to the vector bosons of the theory, V£ . They satisfy 
[T A ,T B ] = if ABC T c , where f ABC are the totally anti- 
symmetric structure constants of the gauge group. Then 
results below are written in terms of the invariants: 

(T A T% J = C(i)5i, (2.4) 



t Since a complex scalar can be written as two real scalars, and 
a Dirac fermion as two Weyl fermions, this entails no loss of 
generality. 



Tt[T a T b ] = I{R)8 AB , (2.5) 

f ACD f BCD = C{G)5 AB^ (2g) 

which define the quadratic Casimir invariant for the rep- 
resentation carrying the index i, the total Dynkin index 
summed over the representation R, and the Casimir in- 
variant of the adjoint representation of the group, respec- 
tively. When the gauge group contains several simple or 
U(l) factors labeled a, b, c, . . . with distinct gauge cou- 
plings g ai the corresponding invariants are written C a {i), 
I a (R), and C a (G). The normalization is such that for 
SU(N), C{G) = N and each fundamental representa- 
tion has C(i) = (N 2 - l)/2N and contributes 1/2 to 
I(R)- For a U(l) gauge group, C(G) = and a represen- 
tation with charge q has C(i) = q 2 and contributes q 2 to 
I(R). The results given below will be presented in terms 
of these group theory invariants for the representations 
carried by the scalar and fermion degrees of freedom. 

The computations in this paper are performed in a 
general gauge with a vector boson propagator obtained 
by covariant gauge fixing in the usual way: 

-i5 AB [^/k 2 + (£ - l)fc"fc7(fc 2 ) 2 ]. (2.7) 

Here £ = for Landau gauge and £ = 1 for Feynman 
gauge and £ = 3 for the Fried- Yennie gauge [31], for a 
vector boson carrying 4- momentum . Infrared diver- 
gences are dealt with by first computing with a finite 
vector boson mass, and later taking the massless vector 
limit. All contributions involving gauge boson loops im- 
plicitly include the corresponding contributions of ghost 
loops. 

For each Feynman diagram, the integrations over in- 
ternal momenta are regulated by continuing to d = 4 — 2e 
dimensions, according to 

J d A k -> (27r^) 2e J d d k. (2.8) 

In the dimensional regularization scheme, the vector 
bosons also have d components, while in the dimensional 
reduction scheme they have d ordinary components and 
2e additional components known as epsilon scalars. For 
the present case of massless vector bosons, this means 
that the 4-dimcnsional metric in the vector propagator 
of eq. (2.7) is replaced by 

rj^/k 2 -> g^/k 2 +g^/(k 2 + m 2 ), (2.9) 

where g^ u is projected onto a formal c?-dimensional sub- 
space, and g^ onto the complementary 2e-dimensional 
subspacc. Counterterms for the one-loop sub-divergences 
and the remaining two-loop divergences are added, ac- 
cording to the rules of minimal subtraction, to give finite 
results, which then depend on the renormalization scale 
Q given by 

Q 2 = 47re-V- (2-10) 

Logarithms of dimensionful quantities are always written 
in terms of 

\nX ee ln(X/Q 2 ). (2.11) 
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The resulting renormalization schemes are known as MS 
[32] and DR [33] , respectively, for the cases in which g**" 
is not and is included. 

The epsilon-scalar squared mass parameter m 2 appear- 
ing in the DR scheme is unphysical. One could set to 2 
equal to zero at any fixed renormalization scale, but then 
it will be non-zero at other renormalization scales, since it 
has a non- homogeneous beta function [34] . Furthermore, 
under renormalization group evolution it will feed into 
the ordinary scalar squared masses in the DR scheme. 
Fortunately, a redefinition (given in [35] at one- loop or- 
der, and at two- loop order in [36]) of the ordinary scalar 
squared masses completely removes the dependence on 
the unphysical epsilon scalar squared mass m 2 from the 
renormalization group equations and the equations re- 
lating tree-level parameters to physical observables. The 
resulting DR scheme [35] is therefore appropriate for 
softly-broken supersymmetric theories such as the Min- 
imal Supersymmetric Standard Model (MSSM). In this 
paper, calculations will be presented simultaneously in 
all three schemes, using the following two devices. First, 



1 for MS 

for DR, DR'. 



(2.12) 



Second, terms that involve the unphysical parameter m\ 
should be construed below to apply only to the DR 
scheme, not the DR or MS schemes. 

A main objective of this paper is to compute the two- 
loop scalar self-energy 



1 n« 



16tt 2 



(16tt 



n^ + ..., (2.13) 



2\ 2 »J 



a (complex, in general) symmetric matrix, as a function 
of 



s = —p 



(2.14) 



where p M is the external momentum. Note that s is 
taken to be real with an infinitesimal positive imaginary 
part to resolve the branch cuts. The self-energy function 
Hij is calculated as the sum of connected, one-particle 
irreducible, two-point Feynman diagrams. It is gauge- 
dependent, but can be used to obtain a gauge- invariant 
physical squared mass, defined as the position of the com- 
plex pole, with non-positive imaginary part, in the prop- 
agator obtained from the perturbative Taylor expansion 
of the self-energy function. For scalar particles with tree- 
level renormalized (running) squared masses m\ , the two- 
loop pole squared masses 



Sk 



M 2 - iT k M k 



(2.15) 



are obtained as the solutions to 

Dot [(to 2 - 8k )8ij + U tJ (s k )] = 0. (2.16) 

A gauge- invariant and renormalization scale invariant so- 
lution at two-loop order is obtained by first expanding 



the self-energy in a Taylor series in s about the tree-level 
squared mass, with the result for the complex pole mass: 



Sk = 



m 2 k + 



1 



(i) 



IQ^kk 



+ 



1 



(16tt 2 ) 2 



n (2 ) -i- tt ±; tt 

ll kk t Hl.l. n 



'-kk 



kk kk 



-5>£>) a /(m2-m?)' 



(2.17) 



where the prime indicates differentiation with respect to 
s, and all self-energy functions on the right-hand side are 
evaluated with s — > m\. This assumes that, as is the case 
for example for sfermions in the MSSM, the scalars that 
mix with each other are not degenerate, so that the last 
term is a well-defined part of a perturbative expansion. If 
(nearly) degenerate scalars do mix, then the appropriate 
version of (nearly) degenerate perturbation theory should 
be used instead. 

One can also obtain a solution iteratively, by first tak- 
ing Sk — m\ as the argument of the self-energy eq. (2.16), 
and then taking the resulting value for Sk and substitut- 
ing it in as the argument of the self-energy function, re- 
peating the process until sufficient numerical convergence 
is obtained. In this case, since Sk has a negative imag- 
inary part on the physical sheet and the argument s of 
the self-energy is taken to be real with a positive imagi- 
nary part, the self-energy for complex s can be defined in 
terms of its Taylor series expansion about a point on the 
real s axis. However, for a theory with massless gauge 
bosons, the terms of a given loop order in the expansion 
of the self-energy have branch cuts. For example, at one 
loop order, 



n«( S ) = (£ - 3)g 2 a C a (z)S t3 m 2 P( S /m 2 ) + ... (2.18) 



where 



P(x) = (x- l/x)ln(l 



l£) 



(2.19) 



and the ellipses refers to terms without branch cuts. This 
yields a result that is not perturbative in the gauge cou- 
pling, unless one takes the Fried- Yennie gauge-fixing con- 
dition £ = 3. So, although the pole mass is formally 
gauge invariant, this iterative procedure has quite poor 
convergence unless £ = 3. At least in the examples given 
below in section V, I find that the iterative procedure in 
Fried- Yennie gauge gives good agreement with cq. (2.17), 
the difference being formally of three-loop order in any 
case, and the implementation of eq. (2.17) is simpler and 
computationally faster. The checks of gauge invariance 
and renormalization scale invariance for particular ex- 
amples are also obtained most straightforwardly from 
eq. (2.17). 

The results below will be written in terms of two- 
loop integral basis functions, following the notation given 
in [37, 38]. The one-loop and two- loop integral func- 
tions are reduced using Tarasov's algorithm [47, 48] to a 
set of basis integrals A(x), B(x,y), I(x,y,z), S(x,y,z), 
T(x,y,z), U(x, y, z, u), and M(x,y,z,u,v), correspond- 
ing to the Feynman diagram topologies shown in fig. 1. 



4 





T(x,y,z) 




S{x,y,z) 




U{x,y,z,u) 




M(x, y, z, u, v) 



Here x, y, z, u, v are squared mass arguments, and the 
arguments s and Q 2 are not shown explicitly, because 
they are the same for all functions in a given equa- 
tion. The name of a particle stands for its squared mass 
when appearing as an argument of a loop-integral func- 
tion. A prime on an argument of one of these functions 
indicates a derivative with respect to that argument, 
so that T(x,y, z) — —S(x',y,z) = —dS(x,y,z)/dx. 
Also, I(x,y,z) = S(x,y, z)\g—o. It is often useful to 
define the functions V(x,y, z,u) = —U(x,y',z,u), and 
T(0,y,z) = ]im x ->o[T(x, y, z) + B(y, z)lnx]. These and 
B(x,y') can be reduced to combinations of other basis 
functions, but they arise quite often in applications in 
such a way that explicit reduction would needlessly com- 
plicate the expressions. The basis integrals contain coun- 
terterms that render them ultraviolet finite. The precise 
definitions, and the calculation of these functions and a 
publicly available computer code for that purpose, are 
described in [37, 38]. 

The one-loop and two-loop Feynman diagrams in the 
approximation used in this paper are shown in fig. 2, 
labeled according to a convention described in detail in 
ref. [26] . The results of their evaluations are described in 
the next section. 



FIG. 1: Feynman diagram topologies for the one- and two- 
loop self-energy basis integrals used in this paper. The dot on 
the T diagram means that T(x,y,z) = — dS(x, y, z)/dx. The 
precise definitions of the integral functions and methods for 
their evaluation are described in [37, 38]. 



III. TWO-LOOP SELF-ENERGY FUNCTIONS 

I begin by reviewing the result at one-loop order. Using 
the notation for loop integral functions used in [26], the 
one-loop self-energy function matrix for real scalars is: 



n (i) = }_ X « kk A s (m 2 ) + ^ kn \> kn B ss ( m lml)+Re{y Km y KNj }B FF (ml,m 2 N ) 

+Re[y Km y K ' N ' j M KK ,M NN ,]B^(m 2 K ,m 2 N ) + 5 ij9 2 a C a (i)[Bsv(ml 0) + A v {0)] (3.1) 

Here I have not written the contribution from massive vector bosons, which can of course be consistently included in 
the one-loop part even if it is neglected in the two- loop part as below. In ref. [26], the results were given in the MS 
and DR schemes, for which Ay(0) vanishes. In the DR scheme, for massless vector bosons, one has instead: 

A y (0) = -2m 2 . (3.2) 

Going from the DR scheme to the DR scheme at one loop order just removes this term [35]. Similarly, at two- loop 
order, the change of scheme given in [36] removes all terms that depend explicitly on the unphysical epsilon-scalar 
mass. The other functions, due to the first four diagrams in fig. 2, are given in terms of the basis functions by: 

A s (x) = A(x) = x(Tnx-l), (3.3) 



1 



-B SS {x, y) = -Bpp{x, y) = B(x, y) = - dt ln[tx + (1 - t)y - <(1 - t)s - is] (3.4) 

1 Jo 

B FF (x,y) = (x + y-s)B{x,y)-A{x)-A{y) (3.5) 

B S y(x,0) = (3-t)(x + s)B(p,x) + (3-2£)A(x)+2(Z-l)s (3.6) 

= is + 3x-3xBx + ^-3){s-x)^x+(l + x/s)bx{l-s/x-ie)]. (3.7) 
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FIG. 2: The one-loop and two-loop Feynman diagrams for scalar boson self-energies in the approximation of this paper. Dashed 
lines stand for scalars, solid lines for fermions, and wavy lines for massless vector bosons. Diagrams involving vector boson loops 
also include the corresponding ghost loop diagrams. The label for each diagram refers to a corresponding function obtained as 
the result of the two-loop integration. All counterterm diagrams for each diagram are included in these functions, rendering 
them ultraviolet finite. For each diagram with a fermion loop, fermion mass insertions (indicated by adding a bar to the 
corresponding subscript F in the name) are to be made in all possible ways. Diagrams indicated by (*) vanish identically in 
the MS scheme, but not in the DR scheme with non-zero epsilon-scalar masses. Ref. [26] explains the naming convention. 
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At two-loop order, one can write the self-energy function as a sum of contributions from diagrams with 0, 1, and 2 
or more vector lines: 



n 



(2,0) 



+ n (2>1) + 

y ^ y 



(2,2) 



(3i 



First, consider 11^ •' , the two-loop self-energy contributions from diagrams with two or more massless vector boson 
(or ghost) propagators. The pertinent individual diagrams shown in fig. 2 contribute with group theory factors 
proportional to: 



Vvsssv,Yvssv ■ 


g 2 a g 2 C a (i)C b (i) 


(3.9) 


ZvSVS, Uyssy, Ssw ■ 


9 lC a (i)[glC b (i)-\glC a {G)] 


(3.10) 


M-vssvs ■ 


glC a (iMC h (i)-\glC a {G)} 


(3.11) 


Vsvvw, Wyyyy , Y S VVV , Xyyy, MyvSSV ■ 


g 4 a C a (i)C a (G) 


(3.12) 


VSVVFF, V SVVTF' WyyFF, ^VVFF : 


9a,C'a{i)Ia{Rfcrmions) 


(3.13) 


VSVVSS, WyysS, ^SVVS,^VVS ■ 


9a C<i (i )Ia (-Rscalars) • 


(3.14) 



Reorganizing the contributions in terms of the four independent group theory factors, and writing them in terms of 
the basis functions, I obtain: 

l4 2 ' 2) = Sii&CaWltiCMF^ (3.15) 

with loop integral functions defined by: 

F x (ar) = -(3a; + s) 2 M (0, x, x, 0, x) - (12.x + 4s)[/(x, 0, x, x) + 8xU(x, 0, 0, 0) + (2s - 6x)T(x, 0, 0) 
-3S(x, x, x) + [(3a; + s)B(0, x) + 4s - 12a; + (6 + 2s/x)A(x)]B(0, x) + 7A(xf/x 
-19A(x) + lOx + 15s/2 + 4[4x - 3A(x)](x/s) ln(l - s/x) 

+ (1 - £){-8xU(x, 0, 0, 0) + 4(s + x)T(x, 0, 0) + 4[xB(0, x) + (1 + s/x)A(x) + 3x + s]B(0, x) 
+A(x) 2 /x - 5A(x) + Tlx + 2s + 2[4x - 3A(x)](x/s) ln(l - s/x)} 

+ (1 - { -2xU(x, 0, 0, 0) + (s + a;)T(a;, 0, 0) + [xB(0, x) + (l + s/x)A(x) + 4x + 2s]S(0, x) 
+A{xf/x + 2A(x) + 5x/2 - 3s/2| + <%g[s - 2x + 4A(ar)] 

+m 2 [2A(x)/x - 2 - (8x/s) ln(l - s/x) + (1 - £){-2A(a;)/a; - 6 - (4x/s) ln(l - s/a;)}] (3.16) 

F 2 (x) = [(3a; + s) 2 /2]M(0, a;, a;, 0, a;) + (s - x) 2 M(0, 0, x, x, 0) + (6x + 2s)?7(x, 0, x, x) + 10x£/(x, 0, 0, 0) 
+ (6s - 10x)T(x, 0, 0) + (3/2)S(x, x, x) + (13s/2 - x/2)B(0, xf 
+ [7s/3 - 77x/3 + 6(1 + s/x)A(x)]B(0, x) + 6,4(x) 2 /x - 47A(x)/3 - lOx + 49s/4 

+ (! - - x) 2 /2]M{0, 0, x, x, 0) + (s/2 - llx/2)T(x, 0, 0) + 5xU(x, 0, 0, 0) 

+ [(3s - x)B(0, x) + (3 + s/x)A{x) - 5s - 29x]B(0, x)/2 + A{x) 2 /x - 7A(x) - 8x + 4s| 

+ (1 - £) 2 {xT(x, 0, 0) - (s/2 + x/2)U{x, 0, 0, 0) + [9x/4 + 5s/4]B(0, x) + 3A(x)/2 + 5x/8 - s/8| 

+<W( S + X ) B (°> x ) + x / 2 ] - 10m2 ( 3 - 17 ) 
F 3 (x,y) = {-4x[{s-xf + 16{s + x)y-24y 2 ]/15y(s-x)}T(x,y,y) 

+{[22s 3 - 14s 2 x - 38sx 2 + 30x 3 - 120s 2 y + 48sxy + 8x 2 y + 128sy 2 - 64xy 2 ]/15(s - x) 2 }T(y, y, x) 
+{[12s 2 x - 4s 3 - 12sx 2 + 4x 3 - 86s 2 y + Usxy + 42x 2 y + 128sy 2 - 64xy 2 ]/15y(s - x) 2 }S(x, y, y) 

+{2(s + x)[A(y)/y + 1] + 2(s - x) 2 /15y}B(0, x) + {[8(s - x) 2 + 32(s + x)y]A(x)A(y) 

+ [4sx - 26s 2 + 22x 2 + 64sy - 32xy]A{y) 2 + [2s 3 - 6s 2 x + 6sx 2 - 2x 3 + 52s 2 y - 40sxy - 12x 2 y 
-32sy 2 + 6Axy 2 ]A(x) + [56s 2 y - 8x(s - x) 2 + 144sxy - 264x 2 y - 256sy 2 + l28xy 2 ]A(y) 
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+5s 4 /2 - 23s 3 a;/2 + 39s 2 x 2 /2 - 29sx 3 /2 + 4x 4 + 20s 3 y - 38s 2 a;y - 32sx 2 y + 50x 3 y 

-244s 2 ?/ 2 + 2AQsxy 2 - \QQx 2 y 2 + 256sy 3 - 128a;y 3 | /15y(s - xf + <%js(-2y) + 2m 2 (3.18) 

F 4 (x,y) = {-x[{s-xf + 56{s + x)y + 96y 2 ]/15y(s-x)}T{x,y,y) 

+{[13s 3 - lls 2 x - 20s 2 y - 17sx 2 - 8sxy - 128sy 2 + 15x 3 + 92x 2 y + 64xy 2 ]/15(s - xf}T(y, y, x) 
+{[(x - sf - 69s 2 y + 2&sxy + 43x 2 y - 128sy 2 + Mxy 2 ]/lhy(s - xf}S(x, y, y) 

+{(s + x)[A{y)/y + 1] + (s - xf/30y}B(0, x) + {[2(s - x) 2 - 32{s + x)y]A{x)A{y) 

+ [32xy - 14s 2 - Usx + 58x 2 - 64sy]A(yf + [(s - xf /2 + 28y(s 2 - x 2 ) + 32y 2 (s - 2x)]A{x) 

+ [4sx 2 - 2s 2 x - 2x 3 + 164s 2 y - WAsxy + Ax 2 y + 256sy 2 - 128xy 2 ]A(y) + 5s 4 /8 - 23s 3 x/8 

+39s 2 x 2 /8 - 29s:r 3 /8 + x 4 + 45s 3 y/2 - 42s 2 xy - 51sx 2 y/2 + A5x 3 y + 24s 2 y 2 - 280sxy 2 

+360x 2 y 2 - 256sy 3 + 128a;y 3 }/15</(s - xf + 6 m 4A(y). (3.19) 

In the limit s — > x, the function Fx has a logarithmic singularity of the form ln(l — s/x), and F3 and F 4 have l/(s — x) 
and l/(s — xf singularities in individual terms. After taking into account the identities mentioned in Appendix A, 
one can check that F2, F 3 , and F 4 are actually finite in that limit. [Because of these same identities, which hold 
between the basis functions due to coincident and vanishing arguments, the representations given in eqs. (3.16)-(3.19) 
are not unique.] Also, the functions F\ and F 2 are both dependent on the gauge-fixing parameter £. It is therefore 
useful to define, for the purpose of computing the finite and gauge-invariant pole mass, the functions: 

Fi(a;) = lim[F 1 (x) + [B sv {x,0)+A v {0)]B' sv (x,0)} 

S— >X 

= xfl27r 2 - 16tt 2 In 2 - ^ + 24((3) - ^rnz + ^-hfx) + 8 m x(4\nx - 5) + m 2 (6\Ex - 4) (3.20) 
V 8 2 2 / 

F 2 (x) = limF 2 (x) 

S—fX 

/1 147 107r 2 40Q 1Q 9 \ / Q \ 

= x(^-^-+87r 2 ln2-12C(3)-^ln 2 ;+^ln x) +6^x1 | - 21narJ -10m 2 (3.21) 

F 3 (x,y) = lim F 3 (x,y) 

S— *X 

2 

= 2{x - y 2 /x)U 2 (l - x/y) + 8(x - y)f{^Jy~Jx) + y(\8 + ^ - 6rna; - 6Ey) 
/ 49 7r 2 19 2 \ 

+ x {-— - y + Y lnx ~ 2lnxln y + ln y) + <W- 2 2/) + 2to ' (3-22) 

Fi(x,y) = lim F 4 (a;,y) 

= (x + 6y + y 2 /x)Li 2 {l - x/y) + 8(x + y)f(y/y/x) + y[-A - tt 2 (1 + y/Qx) + 7ln(x/y) + 3rn" 2 y] 
/ 75 7r 2 25 1 2 \ 

(-— - — + —\nx - \nx\ny + -In yj + ^^(lny - 1) (3.23) 

f(z) = z{Li 2 ([l-z]/[l + z])-Li 2 ([z-l]/[l + z]) + ^ 2 /4}. (3.24) 

I 

The reason for including the term proportional to +^ms( — 2x ) + 2m 2 , (3.25) 

B' sl/ (x,0) = dBsv( x i 0)/ds in the definition of F\ in „ /87r 2 107 25 5 2 \ 

eq. (3.20) is that this includes the appropriate contribu- F 4 (x, x) — x y~^ g — I" ~q^ ux 2^ n x ) 

tion from the U^U^' term in the formula for the pole ^ ^ x (\^ x \\ (3 26) 

mass, eq. (2.17), exhibiting the simultaneous cancellation MS ^ \ ■ ) 
of the gauge dependence and the logarithmic singularity anc j f or y = 0: 
in the limit s — > x . 



+x 



where 



Some useful limits arc, for x = y: F ^ Xj ) = x (yln.T - hA - ^ - + 2m 2 , (3.27) 

~ /23 17 2 \ ~ /25 1 2 75 ir 2 \ 

F 3 (x,x) = x(^— —\nx-lnx) F 4 (x,0) = x(^—\nx - -In x - — —J, (3.28) 
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and for x = 0: 

F 3 (0,y) = 2/(4 - 12my) + ^(-2^) + 2m 2 , (3.29) 
F 4 (0,y) = y(ll + 3bi 2 y) + <^4y(my-l). (3.30) 

Next, consider IT^ 2 ' 1 ), the contributions from diagrams 
with one vector line. These were already given in sections 

I 

ng' X) = \\^ kk gl{c a {k)G s { m l)+C a (i)Gsss{nilml m l)) 

+ \x kn V kn 9l{\\Ca{k)+C a (n) C a {i)}Gss(ml,m 2 n ) 
+[C a (i) - C a (k) + C a (n)]G S sss(ml m 2 , m|, m 2 )} 

+Rc[y XAr W J ]5a{^[C a W + C„(iV) - C (i)]G fF K,m^) 

- C„(A-) + C a (N)]G S SFF(m 2 ,m 2 ,m 2 K ,m 2 N )} 
+Rc{y Km y K ' N 'iM KK ,M NN ,}g 2 a fyc a (K) + C a (N) - C a (i)]G^{m 2 K ,m 2 N ) 



+[C a (i) - C a (K) + C a (N)]G ssln? (ml m lm 2 K ,m 2 N )} (3.31) 

The functions Gss, Gff, Gjrp, Gssss, Gssff, G ss jrp, were defined in equations (5.31)-(5.36) of [26], referring to 
earlier results in that paper. Also I define here: 

G S (y) - Wsssv(y,y,y,0)+X ssv (y,y,0), (3.32) 

= y(-12 + llIny-3ln~ 2 y)-2TO 2 i^/ (3.33) 

Gsss(x,y, z) = Y V sss(0,x,y,z) (3.34) 

= A(z)[(3-2Z)A(x) + (8 + x)(3-Z)B(0,x)]/(x-y) + (x~y), (3.35) 



in terms of functions appearing in eqs. (5.1), (5.2), (5.4), and (5.5) of ref. [26]. All of these functions were written in 
that paper in the MS and DR schemes in a notation consistent with eq. (2.12) of the present paper. To obtain the 
DR scheme results, one should add — 2m 2 lnx to eq. (5.2), and 2m 2 B(x,y') to eq. (5.3), and —2m 2 to eq. (5.17), all in 
ref. [26]. 

The functions Gs{x), Gss{x,y), GFF{x,y), and G-pp(x,y) are each independent of the gauge-fixing parameter £. 
The other functions are are not. However, for the purpose of computing the two-loop pole mass, one can define the 



gauge- invariant functions: 

G S ss{x,x,z) = lim[Gsss(x,x,z)+A s (z)B' sv {x,0)], (3.36) 

Gssss(x,x,y,z) = lim [G S sss{x, x,y, z) + B S s(y,z)B' sv (x,0)}, (3.37) 

G S sFF(x,x,y,z) = lim [G S sFF{x,x,y,z) + B FF (y, z)B' sv (x,0)} , (3.38) 

GsSFF( X > X >y> Z ) = lim [ G SSFF( X > X >y> Z ) + B FT(yi Z ) B 'sv( X '°)} ■ ( 3 - 39 ) 



The first two arguments are taken equal in these functions, because one can consistently neglect the off-diagonal 
entries in the two-loop part of the self-energy when computing the two-loop pole mass. As before, the reason for 
including the terms involving B' sv (x,0) = dBsv{x, 0)/ds is that this naturally includes the corresponding parts of 
the term involving IlWllW in eq. (2.17). 

The result of taking the limits in eqs. (3.36)-(3.39) is: 

G SS s(x,x,z) = {4-3m~x)A(z) (3.40) 
Gssss(x, x, y, z) = 4(x -y + z)M(0, z, x, y, z) + W {x, 0, z, z) + 2U (z, y, z, x) + 4T(0, y, z) - 45(0, y, z)/x 



IV.C, IV.D, IV.E and V of ref. [26], but the results can be 
rewritten in a somewhat nicer form in the case of massless 
vector bosons: 
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-2yT(y, 0, z)/x + (2 - 2z/x)T(z, 0, y) + 2I(x, y, z)/x + 3I(x', y, z) + z(3\nz - 7)B(y, z') 
+2{\Ex + TRz- 8)B(y, z) + (2yEy + 2zTRz - 7x/2 -Ay- Az)/x 



Gssff( x > x > V' z ) 
GsSFF(x,x,y, z) 



= -2G S sss(x, x , y, z) + 6zB(y, z')(l - lnz) + 5 m [-4zB(y, z' 



(3.41) 
(3.42) 



4[(x - yf - z 2 ]M(0, z, y, x, z) + A{x -y- z)U{x, 0, z, z) - AyU{z, y, z, x) 
+4(x -y- z)T(0, y, z) + (x - 2z) (x - y - z)T(z, 0, y)/x + 2y(y + z)T{y, 0, z)/x 
+(4(?/ + z)/x - 1)5(0, y, z) + 2S(x, z, z) + 3{x- y - z)I(x\ y, z) 
+ [1 - 2(y + z)/x]I(x, y, z) + 2(x -y - z)z(3hiz - 5)B(y, z') + [15y + 17 'z - 7x 
—2(x + y + z)hix + (x — y — 3z)lnz]_B(y, z) + 31n.x(ylny — y + zlnz — z) 

-y[3 + 2(y + z)/x]Tay + z[5hiz - 15 - 2(y + z)/x]hiz + 4(y + z f/x - 15.t/4 + lly/2 + 33z/2 
+<%Tsfe _ x ~ z ) B (y> z) + 2(x-y- z)zB(y, z) - yhiy + zhiz - 3x/4 + 2y + 2z\ (3.43) 

Finally, consider the contributions from diagrams without any gauge interactions, n^'°\ These were already given 
in sections IV. A and IV. B of ref. [26] in exactly the same notation, and so will not be repeated here. However, for 
the purpose of computing the pole mass, it is convenient to define: 



ff (2 ' 0) = lim 



n 



(2,0) 



n<;>{i(A< 



kn\2r>l 



.KNi\2t3I 



B FF (m K ,m N ) 



+Rc[ y Km y K ' N ' l M KK ,M NN ,}BL T (m 2 Kl m 2 N )} 



(3.44) 



where the primes indicate differentiation with respect to s. This incorporates the rest of the terms involving n^n^' 
in cq. (2.17). (The derivatives of one- loop functions with respect to s arc easy to obtain analytically, and are given in 
the present notation in ref. [26].) 
The pole mass is now obtained as: 



16tt 2 



n 



(i) 



(16^ 2 ) 2 



n f ,o) + fi ( 2 , 1) + ri ( 2 , 2 ) + E(ri a) )2/(TO? _ m?) 



(3.45) 



with no sum on i implied. Here is obtained 

(2 1) 

from ' by replacing the functions Gsss, Gssff and 

~ ~ ~ ~ (2 2 ) 

Gssff with G sss, Gssff, and G SSFl? , and TL) ' is 

(2 2) 

obtained from ' by replacing the functions i*i,2,3,4 
with Fx.2,3.4, and 

n& } = lim n^. (3.46) 

It is a nice check that the limit s — > x gives a finite 
result for the pole mass here. Also, the independence of 
the pole mass with respect to the choice of gauge-fixing 
parameter, up to terms of three-loop order, now follows 
immediately from the absence of £ in eqs. (3.20)-(3.23) 
and (3.40)-(3.44). Note that this relies on cancellations 
involving the two-loop and iterated one-loop self-energy 
function contributions to the pole mass. 

IV. TWO-LOOP SUSYQCD CORRECTIONS TO 
SQUARK SELF-ENERGIES AND POLE MASSES 

As an example application of the preceding results, 
consider the two-loop strong (SUSYQCD) contributions 
to the squark masses in supersymmetry. Consider an ap- 
proximation in which the squark mixings respect family 



symmetry, but can mix left- and right-handed squarks. 
(This is a slight generalization of the usual assump- 
tion that only the third-family squarks have a significant 
mixing.) The tree- level squared- mass matrices for each 
squark flavor can then each be treated as 2 x 2 matrices 
m? in the gauge eigenstate basis. They are diagonalized 
by unitary transformations: 

(£)-*(*)• , "» 

for q = u,d, c, s, t, b, chosen so that 

Thus, gi, <?2 are mass eigenstates with squared masses 
m Ii' m %' wnne Ql, Qr are the gauge eigenstates. Uni- 
tarity of the matrix X q allows one to write 

where L qi = R*- 2 = Cg, and R qi = -L*- 2 = Sg, with 

Kf + kf = l- (4.4) 
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If the off-diagonal elements of the squared mass matrix 
are real, then Cq and Sq are the cosine and sine of a squark 
mixing angle for each of q — t, b. Also, to a good approx- 
imation in most realistic models, s q — 1 and c q = for 
q = u,d, c, s. For convenience, I define the following com- 
binations: 

Pi = L q R q //,; . (4.5) 

Ni = Lq^m+Rq^ (4.6) 



when q~i and q~j are of the same flavor, and Pf = Nf = 
otherwise. (Here, squarks are complex scalars, so the 
heights of indices are significant.) Also, in the following 
C q = 4/3 and C G = 3 are the quadratic Casimir invari- 
ants of the fundamental and adjoint representations, and 
I q = 1/2 is the Dynkin index of the fundamental. 

The self-energies will then likewise be 2 x 2 matrices 
for each of the 6 flavors. The SUSYQCD contribution to 
the one-loop self-energy in the DR scheme is: 



nf b = glC q P*Plk s {q k ) + 2S>B FF (q, g) - 2N 1 i m q m s B TF (q, g) + <5?'B SV 0) 



(4.7) 



where k is summed over the two mass eigenstates of the same flavor as (The other, non-SUSYQCD, one-loop 
contributions can be found in ref. [30].) 

The contributions from two-loop SUSYQCD diagrams with no gluon propagators are: 

Of = AgjC 2 \-P^P^PlP^sss{q k An, q m ) + \p* SsFF^Qk, 9, §) 



1 

2' 

-5 3 - Ntm„m 



^ P i N r p i m q m g W sSFF^k,q m ,q,9) + {Lq z L* qj \Lq k \ 2 + Rq t R qj \Rq k \ 2 )V FFFFS (g,q,q,g,qk) 

i^k " l q" l gV FF TFs(9, Q, Q, 9, Qk) + { L qi R Sj R qk L *qk + R V L *q 3 L qk R *q k ) m g^F F FFs(9 ' 9> 1, 9, Qk) 

+{ L m L *q j \ R Sk ? + R qi R *qi \ L q k 1 2 ) m \ ^fTffs <?> Q, 9, Qk) - N? m q mgV F FF F s (g , q, q, g, q k ) 

+ ( L q, R *q ] L q k R q k + R q* L % R q k L \ k ) m \ m \ V FFFFS (ff> <?> 9' 9, Qk) 

+gi[AC 2 q - 2C G C q ] [{L qi Rt.R qk L* qk + Rq i L*q. Lq k )M FFFF g (q, g, g, q, q k ) 
-N?m q m s M FFFF - s (q, g, g, q, q k ) + {L qi L*-. \L qk \ 2 + R qi R~. \R qk | 2 )™|M fWfs ( (? , g, g, q, q k ) 
-5jN^m q m g M FFFFS (q,g,g,q,q k ) + (L qi L* q . \Rq k \ 2 + R qi R* q . \L qk \ 2 )m 2 M FF -p FS (g, q, q, g, q~k) 

+(L qi Rt.L qk Rt k +R qi LlR qk Ll)m 2 q mlM F ^ FFS (q,g,~9,q,qk) + -^P^P^P^SSsiQkAmAn) 
+AglC q Iq SfV FFFFS (q,g,g,q r ,qr) + Sfm 2 gV FFTFS (q,g,g,q r ,qr) - 2N J i m qr mgV FFFFS (q,g,g,q r ,q r 



+Nf N^m q m qr V FFFFS (q, g, g, q r , q r ) - 26jN^m qr m g V FF - FFS (q,g,g,q r ,q r ) 
+N^N^m q mq r m 2 g V FFFFS (q, g, g, q r ,q r ) 



-P?P J k P:P:Ssss(qk,qr,q s ) 



For two-loop diagrams with one gluon propagator, 

nf = 9 tSf{CqC G [G FF (qrg) + GssFF(<ji, q%, q, g)] + [2C q - C G }C q G S sFF(qi,qi,g,q)} 

-glN 3 i m q m g {C q C G [G TF {q,g) + G SSFF (qi,qj,q,g)] + [2C q - C G ]CgG SSFF (<z;, q h g, q)} 
+9tC 2 q P*Pl [G s (q k ) + GsssiQuQj^k)} ■ 
Finally, for two-loop diagrams with two or more gluon lines, 



n 



(2,2).? 



33 



sic, 



i w 9 



CqF^qt) + C G F 2 (qi) + C G F 3 (q l7 g) + I q ^{P 3 (ft, Qr) + ^(ft, 9r)} 



(4- 



(4.9) 



(4.10) 



In eqs. (4.8)-(4.10), indices r,s are used when all 12 left-handed and right-handed quarks and 12 squark mass eigen- 
states should be summed over. Indices k, m, n are used when only the two squarks with the same flavor as the external 
particle states q~i, q~j are summed over. The quark corresponding to the external squark flavor is called q. 

The squark pole squared mass M|. can now be obtained following the discussion of the previous section; see 
eq. (3.45). This yields: 



Mi 



to? + -—Il (1)t + 



(16^ 2 ) 2 



nf.o) + nf^) + np + ^npnWVK-m; 



(4.11) 
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with no sum on i implied. Within the approximation specified above for sfermions in the MSSM, j can actually take 
on only one value for a given i, and the sfermions that mix are always non-degenerate. Here HQ is obtained from 



ll- 1 ^ by taking s 

by Gsss, Gssff, and G SSFF , and IT? 2 ' 2 ^ is obtained from H^ 2 ' 2 ^ by replacing the functions ^1,2,3,4 by ^1,2,3,4, and 



Also, II. 2 ' 1 ^ is obtained from II- 2 ' 1 ^ by replacing the functions Gsss, Gssff and G SSFF 



n 



(2,0) 



lim 



% 2 C q {B' FF (ir9) 



(4.12) 



I have checked that this result for the squark pole mass is invariant under rcnormalization group evolution of the 
parameters, up to terms of three-loop order. That check is somewhat messy when non-trivial mixing is involved, and 
so will not be presented here in the general case. Instead, it will be shown explicitly in two simplified special cases in 
the next section. 



V. SIMPLE EXAMPLES 

A. Squarks without mixing 

In this subsection, I consider a simple (probably non-realistic) example, to demonstrate the typical size of the 
two-loop contribution to the squark pole masses. As an approximation, suppose that all squarks are degenerate in 
mass, with therefore no mixing, and that all of the quark masses can be neglected. Then, reorganizing the results of 
the previous section by common group theory factor rather than number of gluon propagators, one obtains for the 
pole squared mass of a given squark qf. 



.^L_ S |C, [C q n^ (q u g) + C G TlW {q u ~g) +I q Y^ lK 2c >(ft, g, q r )] +... (5.1) 



where 



U^(x,y) = As(x) + 2B FF (0,y)+B sv (x,0) (5.2) 
n (2a) {x, y) = X sss (x, x, x) + S sss (x, x, x) + 4yM FFFFS (0, y, y, 0, x) + AY FFFFS {y, 0, 0, y, x) 
+2MV S sFF(x,x,0,y) + G s (x) +G S ss(x 1 x 1 x) + 2G SS ff(x, x, y, 0) + Fi{x) 

+2Il«(.x,y)B' FF (0,y) (5.3) 

n {2b) (x,y) = --S SS s{x,x,x) - 2yM FFFFS (0,y,y,0,x) + G FF (0,y) + G S SFF(x,x,0,y) 

-G SS ff(x, x, y, 0) + F 2 (x) + F 3 (x, y) (5.4) 

& 2c) (x,y,z) = S S ss(x,z,z) + 4V FFFF s(0,y,y,0,z) + 4yV FFFFS (Q,y,y,0,z) + F 3 (x,0)+F 4 (x,z). (5.5) 

All of the basis integral functions on the right-hand sides of these equations are to be evaluated at s = x. This 

example has the virtue that all of these functions can be given analytically in terms of polylogarithms, with the result 
(in the DR scheme): 

n (1) (x,y) = 2x[l +1n(y/x) + (1 - y/x) 2 ln(l - x/y)] + 6y - Ayhiy (5.6) 
& 2a 1 (x, y) = 8(x - y) 2 M(0, 0, x, y, 0) - 8(x - y)yM(0, y, y, 0, x) + (24a; - 8y - 12y 2 /x)U 2 (l - x/y) 
+ (1 - y/x) 2 (2x + 4y- 4y 2 /x) ln 2 (l - x/y) + 4(1 - y/x)[6x - 2y + 4y 2 /x - (x + y)Tnx 
+(x -y- 2y 2 /x)lKy] ln(l - x/y) + (Ux - Ay) ln(x/y) + 8[y]n(x/y) + 3x - y + y 2 /x}hiy 

+24(y - x)Tnx + [24((3) - 6 + 16tt 2 (1 - In 2)]x - (60 + 87r 2 /3)y + (2tt 2 - I2)y 2 /x (5.7) 
fi< 26 > (x, y) = A{x - y) [{x + y)M(0, x, y, 0, y) + yM(0, y,y,0, x)} - 2(x - y) 2 [2M (0, 0,x,y,0) + M(0, 0, y, y, 0)] 
+ (2x - I2y + \2y 2 /x)U 2 {\ - x/y) + 8{y - x)[f(y^) + (1 - y/x) ln 2 (l - x/y)] 
+ [(10x - 8y)TEx + (44y - 16a; - 30y 2 /x)Ey + 19a; - 80y + 6ly 2 /x] ln(l - x/y) 
+ (2y — x)ln x + (8x — 4y)lna;lm/ + (20y — 7x)hi y — (25x + 16y)lna; + (19a; — 66y)\ny 



12 



+x[21 - 12C(3) - 26tt 2 /3 + 8tt 2 In 2] + y(123 + 14tt 2 /3 - 2n 2 y/x) (5.8) 
n (2c) (^?A z) = [(x - y)(y - z)(xy - by 2 + 3xz + yz)/xy 2 ]{2Li 2 ([y - z]/[y - x}) + ln 2 (l - x/y)} 
+ [2(y - z)(5y - z)/x]U 2 (l - y/z) + [8z - 4z(x + z)/y + 6xz 2 /y 2 ]U 2 (l - x/z) 
+8(x + z)f{^/zjx) + [6z + 2yz/x - I6y 2 /x + (4z + 12yz/x - 2z 2 /x - 6z 2 /y) \n{z/y) 
+(10y 2 /x - 2y)my](l - x/y) ln(l - x/y) + {[2z(2x/y + 2z/y - 3xz/y 2 )} \n(z/x) 
+ (1 - z/x)(lOy -7x-7z + 6xz/y)} ln(l - x/z) + (16y + 4z - 6xz/y) \n(y/z) + (x - 6y)\ny 

-x\nx + [4z + (y - z)(5y - z)/x + 2z(x + z)/y- 3xz 2 /y 2 } ln 2 (y/z) + (2x + 26y + 4z)hiz 

+ {9x + 8z) ln(x/z) -7x- 40y - 3z + (2tt 2 /3)(2z(x + z)/y - x - 2z - 3xz 2 /y 2 ). (5.9) 

Ref. [38] gave analytic formulas for the master integral cases M(0, 0, x, y, 0), M(0, y, y, 0, x)\ s = x . The integral 
M(0, x, y, 0, y)\ s =x can be reduced using recurrence relations to results found in [14], [15], and was given in the 
present notation in [38]. Also, M(0, 0, y, y, 0) was originally found in [40] and listed in the present notation in ref. [37]. 
The function f(z) is defined in cq. (3.24) of the present paper. 

The renormalization group scale independence of the pole mass M qi can now be checked. The requirement 
dM?/dQ = amounts to: 



,(i) 



"A 



(2) 



9 2 3 C q Q— n«fe,5) 

g\C q Q^- \c q & 2a *> (q t ,g) + C G & 2b > (q t ,g) + I q £ H (2c) {*,§, q r ) 



+gic q 



dQ 
2 

93 



(i) 



9m| 



r ]n«(*,3) 



(5.10) 



(5.11) 



(with no sum on i), where 
.dX „ 1 



Q 



dQ 



fa 



16tt 2 



(16tt 2 ) 



(2) 



(5.12) 



arc the beta functions of running parameters X = 
gs, m,g, m~.. Equations (5.10)-(5.11) can be checked using 
the results 



8^ 

8 (1) 
B {2) 



gl(-3C G + 2N f I q ), (5.13) 

,9 3 5 [-6C 2 + (4C G + 8C q )N f I q ], (5.14) 

g 2 (-6C G + 4N f I q )m g , (5.15) 

gi\-2AC 2 G + (16C G + 32C q )N f I q }m g , (5.16) 

-8 5 2 C,m? (5.17) 



g\C q [(-80C G + 48C 9 + 48^//,)™? 



Hr 



(5.18) 



with 27V/ the number of quark/squark chiral supcrfields 
(12 in the MSSM), and 



8y, 



(5.19) 



2- 2y/x + 21n(y/x) 
+2(l-y 2 /a; 2 )ln(l-.T/y), (5.20) 



|-n«(x,y) = 4 + 4(y/z-l)ln(l- a; /y) 



d 



-41m/, (5.21) 
= — 32y — 16y 2 /x + I6y \n(y/x) 
+ I6y(l-y 2 /x 2 )ln(l-x/y), (5.22) 



Q— n( 2fc )(a;,y) = 12a: + 164y + 12xln{y/x) 

+(12x - 72y + 60y 2 /x) ln(l - x/y) 
-72yhiy, (5.23) 

-Ax - 52y - 8z + Ax ln(x/y) 

+ (24y - Ax - 20y 2 /x) ln(l - x/y) 
+24ymt/, (5.24) 

which in turn follow directly from eqs. (5.6)-(5.9), noting 
that the master integral basis function M(x, y, z, u, v) has 
no explicit dependence on Q. 

This scale independence thus holds up to terms of 
three-loop order. To illustrate it in practice, consider the 
even more special case that the gluino and all squarks 
have equal running masses at an input renormalization 
scale Qo given by the same value, so that Qo = m g (Q n ) = 
niqiQo)- Figure 3 then shows the scale dependence of 
the squark pole mass as calculated from eqs. (5.1) and 
(5.6)-(5.9). To make this graph, the running parameters 
<?3, m,g, and m q are each run from the input scale Qo 
to a new scale Q, using their two-loop renormalization 
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FIG. 3: Dependence of the calculated pole mass for degen- 
erate squarks on the choice of renormalization scale. The 
input parameters are defined by running gluino and squark 
masses taken to be equal at an input renormalization scale 
Qo = m,q(Q ) = mg(Qo), with qs(Qo) = 0.095. These pa- 
rameters are then run to new scales Q, where the pole mass is 
recomputed. The vertical axis is the fractional change in the 
squark pole mass compared to the squark running mass at the 
input scale. Six squark families are assumed as in the MSSM, 
but quark masses, squark mixing, and non-SUSYQCD effects 
are neglected in this graph. 



group equations (5.13)-(5.18). Here I have put in the 
MSSM values, namely C G = 3, C q = 4/3, I q = 1/2, 
and N f = 6, and taken a s (Qo) = fli(Qo)/47r = 0.095. 
At the scale Q, the pole mass is recomputed, and the 
quantity Mg/mq(Qo) — 1 is shown; in the ideal case of 
an exact calculation the resulting line would be exactly 
horizontal. The two-loop result has a slightly improved 
scale dependence, as expected, but the difference between 
the two-loop result and the one-loop result is actually 
much larger than the scale dependence of the latter. This 
demonstrates that the scale dependence does not give an 
adequate estimate of the theoretical error of the calcula- 
tion. 

The result at Q = Qo is: 



Mi = 



+ 



os/32' 
4tt V 3 
664tt 2 327r 2 ln2 



27 



9 



UtJ V 3 

16C(3) \ 
3 ) 



0.849a s + 1.89a|] 



(5.25) 
(5.26) 



There are no logarithms here, since there is only one mass 
scale, so the result gives some idea of the typical intrinsic 
size of the two-loop corrections. The one-loop correction 
is an increase of order 4% in the pole mass compared to 
the running mass evaluated at itself, while the two-loop 
correction adds an additional amount of order 1%. 

More generally, fig. 4 shows the one-loop and two-loop 
corrections to the squark masses, calculated as above, but 



FIG. 4: The SUSYQCD corrections to the squark pole mass, 
as in figure 3, but now as a function of the varying ratio of 
the tree-level running gluino mass Ms(Qo) = m§(Qo) to the 
running squark masses Qo = mq(Qo)- The squark pole mass 
is computed at the renormalization scale Qo- 



now varying the running gluino mass M 3 (Qo) = rrig(Qo) 
at the fixed renormalization scale Qo = rriq(Qo)- The 
two-loop part of the correction is seen to be largest when 
the squark masses are slightly less than the gluino mass, 
and does not exceed 1% over the indicated range. 

In nearly all realistic models of supersymmetry break- 
ing, M 3 /m q < 1.5 at the TeV scale for the squarks of 
the first two families. For M^(Qo) much larger than 
2mg(Q ), there are large negative loop corrections to the 
squark pole mass from gluino loops. This can happen 
for top and bottom squarks in the MSSM, in which case 
the top and bottom Yukawa couplings and scalar cubic 
couplings must be included to give a reliable result. 



B. Squarks in the supersymmetric limit 

Let us next consider the supersymmetric limit for 
squarks for an SU(n) supersymmetric Yang-Mills the- 
ory with gauge coupling g. The gaugino mass vanishes, 
and Nf flavors of quarks and squarks obtain their masses 
solely from a superpotential 



N f _ 



(5.27) 



Here, Qi are chiral superfields transforming in the fun- 
damental representation, and Q i in the anti-fundamental 
representation of the gauge group. Specializing the re- 
sults of section IV to this case, I obtain two-loop squark 
pole squared masses: 



Mf 



1 



167T 



;H« + 



1 



(16tt 2 ) 



rn 



(2) 



(5.28) 
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where 



n 



(i) 



g 2 C q m 2 (8-4lnm 2 ) 



(5.29) 



nf } = g 4 C q (C q mj [40tt 2 /3 - 16tt 2 In 2 
+24C(3) - 28 - 8mm 2 + 8m" 2 m 2 



+C G mi 



66 - 47r 2 +87r 2 ln2-12C(3) 



-361nm, 2 + 61n m- 



AT/ 



in which 

h(x,y) 



h(x,0) 
h(0,y) 



+2/ 9 ^/i(m 2 ,m 2 ) 



4(x + y)[U 2 (l-x/y)-ir 2 /6] 

+a;[21n 2 y — 4\nx\ny + 121n.x 

+16/(V»7^)-22], 

z(-22 - 4tt 2 /3 + 12ma; - 2hi 2 cc) 
0, 



(5.30) 



(5.31) 

(5.32) 
(5.33) 



with f(z) defined by eq. (3.24). 

The renormalization group invariance of this pole mass 
result now follows from 

/3« = g 3 (~'3C G + 2N f I q ), (5.34) 
= -4g 2 C q m t , (5.35) 
0%] = 9^0,(80, -12C G + 8N f I q )m t . (5.36) 

Since supersymmetry is unbroken in this example, the 
result of eqs. (5.28)-(5.30) must be equally valid for the 
quark pole squared masses as for the squark pole squared 
masses derived directly here. 

C. Gauge mediation of supersymmetry breaking 

As another example application, I show how to re- 
produce the result of gauge mediation of supersymme- 
try breaking to MSSM scalars by specializing the results 
above. In this case, the two-loop order result gives the 
leading effect, found in ref. [49]. (A derivation in terms 
of individual diagrams in Feynman gauge, perhaps useful 
for comparison with the treatment here, was later given 
in [50].) Suppose that there exists a new, vector-like, 
heavy "messenger" quark/squark pair (not part of the 
MSSM). The heavy quark is taken to have a Dirac mass 
toq, and its scalar superpartners have a squared-mass 
matrix of the form: 



7l 2 Q A 

A m 2 Q 



(5.37) 



where A is a supersymmetry-breaking effect. Diagonal- 
izing this mass matrix according to eq. (4.2) leads to 



eigenvalues = m Q ± A, with a mixing angle of 

7r/4. This induces masses for the ordinary squarks of 
the MSSM, which can be treated as massless in leading 
order. From eqs. (4.8)-(4.12), the result for the ordinary 
squark masses is: 



Mi 



4V 



FFFFS 



(o,o,o,g,Q H 



(16tt 2 ) 2 

+4V FFFF5 (0, 0, 0, Q, Q_) + 2S SS s(0, Q+, Q-) 
+2F 3 (0, Q) + F 4 (0, Q+) + F 4 (0, Q_)] . (5.38) 

Then, one can use eqs. (3.29)-(3.30) and the results valid 
for s = 0: 



V F FFFs(0,0,0,x,y) = 



1 



-2yl(0,x,y) 



x-y 

—2xyh\x\ny + 2y(x + y)\ny 
+2x(x + y)lnx — 3x 2 — 3y 2 — 4xy 
S SS s(0,x,y) = -1(0, x, y), 

where 



(5.39) 
(5.40) 



l(0,x,y)_=_(x - y)[U 2 (l - x/y) + (\nyf /2] 

— xlnadny + 2a;lna; + 2y\xvy — 5(x + y)/2. (5.41) 

The result for M|. is equivalent to the one originally 
given in ref. [49]. It is not hard to generalize this to 
the SU(2)l and U(l)y gauge groups to obtain the full 
set of predictions for MSSM squark and slepton masses 
in gauge-mediated supersymmetry breaking. 



VI. OUTLOOK 

In this paper, I have found the two-loop contributions 
to scalar boson self-energies, and thus pole masses, in 
a general gauge theory with massless (or light) gauge 
bosons. These results should apply directly to heavy 
scalars in perturbative models of physics beyond the 
Standard Model, provided the W and Z masses can be 
neglected compared to the dominant mass scales in the 
problem. This is quite likely to be a good approxima- 
tion, for example, for the squarks and sleptons in super- 
symmetric theories, where the difference between the full 
two-loop result and the one reported here is suppressed 
by m 2 z I '™ 2 f crmion multiplied by an expansion coefficient 
that is typically a fraction of unity, as well as a weak 
interaction two-loop factor. 

In section IV, I have given the SUSYQCD contribu- 
tions for squark masses. However, it should be empha- 
sized that the computations of all two-loop contributions 
to all of the the sfcrmion pole masses in this approx- 
imation have been reduced to an exercise (admittedly 
tedious, but certainly amenable to automation by a sym- 
bolic manipulation program) in substitution of running 
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couplings constants and tree-level masses into the formu- 
las here and in ref . [26] , followed by numerical computa- 
tion of basis integrals using a program such as [38]. In 
doing so, the Higgs scalar and the electroweak vector bo- 
son sectors can be treated in an approximation where the 
effects of electroweak symmetry breaking are consistently 
neglected in the two-loop parts. The one-loop part can 
of course be treated exactly using the formulas in [30] . 

A convincing guess as to the likely size of the remaining 
theoretical errors is difficult to obtain. The results of the 
example in subsection VA may suggest that three-loop 
effects on squark masses are usually less than a few tenths 
of a percent, but the limited available data on the conver- 
gence of the perturbative expansion here is not clearly in 
support of this conjecture. Also, the scale dependence of 
the result, although quite mild, has often been seen to un- 
derestimate the theoretical error. It should be noted that 
there will also be substantial sources of irreducible ex- 
perimental error, notably uncertainties in as, the gluino 
mass, and the other superpartner masses. It seems likely 
that global fits to many different observables will be nec- 
essary in order to extract the parameters of the underly- 
ing Lagrangian. Clearly, there will be many challenges to 
overcome to go from future experimental data to a clear 
and precise understanding of the origin of supersymme- 
try breaking. 



Appendix 

In this Appendix, I note the existence of some identities 
for two-loop self-energy integral basis function that are 
useful for deriving some of the formulas of section III. 
These include equations (A.14)-(A.20) of ref. [26], and 

= (s -x)[T(x, 0,0) - U(x, 0,0,0)] + sB(0,x) 2 
+ [s - 3x + (1 + s/x)A(x)]B(0, x) 
-x + 2s-2A(x)+A(x) 2 /x. (A.l) 

Also needed are the values [12, 40] in the threshold limit 

s — > x: 



B(0,x) 
T(x,0,0) 

U(x, 0, x, x) 



2-lnx (A.2) 

tt 2 /3- 1/2 - rnx + (hix) 2 /2 (A.3) 

11 271-2 ni" V 2 / A ,s 

31nx + -In x (A.4) 

2 3 2 v ' 



11 7T 2 — 1—2 

U(x, 0,0,0) = — + — - 31nx + -In x (A.5) 
M(0,x,x,0,x) = [tt 2 In 2 - 3<C(3)/2]/x, (A.6) 

and the pseudo-threshold expansions [43] : 

S(x,y,y) = S (x,y,y) + (l-s/x)Si(x,y,y) 

+ (l-s/x) 2 S 2 {x,y,y) + ..., (A.7) 



where 
S (x,y,y) 

Si(x,y,y) 



-[(x-y) 2 /x}|Li 2 (l-y/x)+7r 2 /6] 
+x(31nx - In x - 3/4)/2 



r^2 



+2/ [In x — 4 + 51ny — lnx — 21nxhu/ 
-(y/2x) ln 2 (y/x)] 
y(l - y/x)[U 2 (l - y/x) + tt 2 /6 
+ ln 2 (y/x)/2] - 5x/8 + (x/2)hix 
+y + y\n(y/x) 



(A.8) 



(A.9) 



S 2 (x,y,y) = Si(x,y,y) 



7 



fiVv/z) - lM . (A.10) 



with f(z) defined in eq. (3.24), and 

So(x,x,x) = x 
S\{x,x,x) = x 



35 11— 3—2 ' 

Y + Y lnx '2 lnx 



3 1- 
8 + 2 1UX 



S 2 (x,x,x) = x [5/4 - tt 2 /8] , 



and 



T(x,y,y) = -—S(x,y,y), 
1 d 

T(y,y,x) = --—S(x,y,y), 
1 d 

T(x,x,x) = --—S(x,x,x). 



(A.11) 

(A.12) 
(A.13) 



(A.14) 
(A.15) 
(A.16) 
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